Testing Gravity with the CFHTLS-Wide Cosmic Shear Survey and SDSS LRGs 
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General relativity as one the pillar of modern cosmology has to be thoroughly tested if we want to 
achieve an accurate cosmology. We present the results from such a test on cosmological scales using 
cosmic shear and galaxy clustering measurements. We parametrize potential deviation from general 
relativity as a modification to the cosmological Poisson equation. We consider two models relevant 
either for some linearized theory of massive gravity or for the physics of extra-dimensions. We use 
the latest observations from the CFHTLS-Wide survey and the SDSS survey to set our constraints. 
We do not find any deviation from general relativity on scales between 0.04 and 10 Mpc. We derive 
constraints on the graviton mass in a restricted class of model. 



I. INTRODUCTION 

Observational cosmology has established the fiat 
ACDM model as the standard model of modern cosmol- 
ogy [II, i, S i, H, B 0, Hi • Within this model the energy 
budget of the Universe roughly goes as follows : 4% comes 
from baryons, 20% from Cold Dark Matter (CDM) and 
76% is in the form of a cosmological constant. Those 
numbers are now known to a few percent but despite 
those tight constraints, the exact nature of CDM is still 
unknown and the high value of the cosmological constant 
-or more generally Dark Energy (DE)- keeps challeng- 
ing our deepest understanding of fundamental physics 
[Oj. The unknown nature of the required last two com- 
ponents might cast some doubts on the foundation of 
this model : the linear cosmological perturbation theory 
in General Relativity (GR). Although GR passes direct 
tests probing the Solar System scales (10-^^m) down to 
the laboratory scales (10~'^m) extrapolating its va- 
lidity on more than 10 orders of magnitude to address 
cosmological scales (lO^^m) is questionnable. GR thus 
has to be checked on all relevant cosmological scales in 
as many ways as possible llj. We present in this paper 
such a test on megaparsec scales using two probes of the 
large scale structures of the Universe : galaxy clustering 
and weak gravitational lensing. 

Our approach is phenomenological and we do not aim 
at developing or at testing specific fully consistent al- 
ternatives to GR. Several classes of theory have been 
studied in details in the literature with the aim at ac- 
counting for either the acceleration of the Universe or 
for Dark Matter or both. The first class includes for 
example, a five-dimensional alternative [l^, [l3| (DGP) 
or the addition of non-linear terms in the Ricci scalar 
to the gravitational action [H, [H, [H, El. The 
second class (for now) includes tensor-vector-scalar the- 



ory [19|, 120, l2l|, 121 l2J, l25|, |26l as a covariant the- 
ory encompassing the MOND theories [13] ■ Finally, the 
last class includes for example the ghost condensate the- 
ory [ll. Generic models that verify Birkhoff's theorem 
were also investigated in [2§| . Instead of exploring these 
new theories, we rather propose simple functional forms 
parametrizing the deviation from gravity on Mpc scales 
and constrain them using current observational data at 
low redshift. If any significant deviation from gravity 
was observed, we would then make the connection with 
theoretical models more explicit. Our key assumption is 
that the deviation from GR we explore can be written as 
a modification to the Poisson equation. As such we are 
testing the Poisson equation on cosmological scales. 

Several authors already tackled the tasks of constrain- 
ing alternate theories of g ravity with available or planned 
observations 130, 31, 
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ular, within the context of linearized relativity, Edery 
[5^ showed that it is not possible to build a relativistic 
theory that reproduces galactic phenomenology (fiat ro- 
tation curves, gravitational lensing, or dynamical mass 
measurements via the virial theorem) without introduc- 
ing dark matter (their key point is that the light defiec- 
tion cannot be made to have the right si gn). Following 
a similar approach, Zhytnikov and Nester [53| consider a 
general metric theory that would reproduce the dynam- 
ics of galactic systems without dark matter and reach the 
same conclusion. They advocate in particular the addi- 
tion of a Yukawa like potential to the usual Newtonian 
one as a generic extension for a metric theory of grav- 
ity. White and Kochanek [s^l then considered the same 
extension to GR gravity, but focussed on scales relevant 
to weak-lensing (note that they only discussed the ef- 
fect on the deflection angle and not on the growth rate 
of structure as we will do below). As they discuss, a 
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generic feature of such hnearized extension to GR is to 
modify the standard Poisson equation relating the den- 
sity field to the grav itational potential. In turn, Uzan 
and Bernardeau 55] considered an alternative modifi- 
cation to the Poisson equation that encompasses DGP 
theories. Shirata et al. [IBl and Sealfon et al. [13] con- 
strained the Yukawa type extension to GR using galaxy 
clustering measurements. Note that Sealfon et al. [53| 
also introduced a power-law extension to gravity. Such 
modifications to the Poisson equation were later studied 
using numerical simulations [i^, |43| • 

Our work is an extension to the latter works [i^, 113, 
[S^ . Issl . [sgI . \Et\ . In particular, we use the latest weak 
gravitational lensing and clustering data to test gravity 
on Mpc scales at low redshift. The two alternate the- 
ories of gravity we constrain are introduced in Sec. [TTl 
before detailing their phenomenology. We then describe 
our data-sets and methodology in Sec. IIIII before present- 
ing and discussing our results in Sec. IIVI 



II. 



MODIFIED THEORIES OF GRAVITY AND 
CONSTRAINT METHODS 

A. Theoretical motivations 



To allow for deviations from general relativity opens 
up many possibilities and a wide range of gravity the- 
ories that we cannot explore exhaustively. We restrict 
ourselves to a few somewhat phenomenological models 
that offer both a physical motivation and observationaly 
tractable cosmological implications. 

We will consider two different models that we intro- 
duce below: the Yukawa and the Uzan-Bernardeau type 
models. 

As White and Kochanek [54] , we follow Zhytnikov and 
Nester [HI] who within well defined postulates presents 
some general arguments regarding the shape of a lin- 
earized metric theory of gravitation. These postu- 
lates comprise post-Newtonian slow motion extension 
((w/c)^ <C 1) and weak gravitational field regime relevant 
to the scale we consider in this work. Their description 
includes forces mediated by massless or massive scalar 
and tensor modes. The metric reduces to 



(1) 
(2) 



5oo = (-l + 2<i>) 

where the potential $ is given by 

$(r) = (1 - a)$(r,0) +a$(r,TO) (3) 



$(r,m) = G y ■ 



r — r' 



Thus, to the usual Newtonian potential ^(r, 0) a Yukawa 
type potential $(r, m) is added corresponding to propa- 
gating massive modes. 

Considering the evolution of over-density, 5{r, t) or its 
Fourier transform, i5(k, the standard linear theory of 



perturbations leads to the Poisson equation relating the 
gravitational potential to the over-density. The Pois- 
son equation writes in comoving coordinates and Fourier 
space [11] 



$(k,a) 



3 H'^Q.r. 



/(k,a) 



(5) 



where a is the scale factor, Hq is the Hubble constant, 
ilmo the present matter density. / is a deterministic func- 
tion equal to 1 for the standard gravity. It can easily be 
shown that the inclusion of a Yukawa type potential as 
in Eq. H leads to 

^ /(k) = 1 - a i— ^. (6) 

1 + (—) 

As such, exploring alternative theories of gravity will be 
equivalent for us to test the Poisson equation. 

Observational constraints on similar models (but with 
slightly different notations) coming from galaxy surveys 
[56l.l57|. as well as using also their non- linear evolution us- 
ing N-body simulations [4^ |43| , have already been stud- 
ied. Note that our notation matches that of White and 
Kochanek [s^l and corresponds to that of Shirata et al. 
[5^ and Sealfon et al. [13] with (a — > ~a, A — 1/to) 
and [4^ with (a —a, A — 1/rs). 

Another set of models captures, in the context of su- 
perstring theories, some brane induced phenomenology. 
In such a scenario, a generic feature seems to be the exis- 
tence of two scales below/above which standard gravity 
is altered [5i,[63,[6ll[62]. The smaller scale (of order a 
millimeter or less), irrelevant to our measurements, cor- 
responds to the existence of Kaluza-Klein gravitons. On 
the other hand, above the branes separation scale, we 
also expect gravity to be altered. Since this scale is ex- 
ponentially larger than the previous one, it becomes cos- 
mological. The use of large scale structures as a probe 
of this form of deviation from standard gravity as been 
advocated by [H, [H, [63l. w ho include the model of [13] . 
Uzan and Bernardeau [55|] focus in particular on weak 
gravitational lensing, as such we will follow their work 
more closely. They describe in real space the violation of 
Newton's law above a given physical scale, mul- 
tiplicative function to the standard Newtonian potential, 
$(r,0), 

<l'(r) = cl.(r,0)^^ (7) 

so that on small scales, r ^ r^, we recover the usual 
Newtonian gravity. As before, this translates into a mod- 
ification to the Poisson equation with a multiplicative 
function, /^^ [13, [11 

/^^(fc,a,r,) = /(fc,a) 



(8) 
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is an arbitrary scale and is the parameter we will con- 
strain later on. 

From now on, we will focus on these two models 
(Yukawa and UB) and study their phenomenological im- 
plications on cosmological scales (~ 10 Mpc) before con- 
straining their free parameters. We deliberately restrict 
our conclusions from those models to the scales we are 
probing, i.e. we do not explore the parameter space where 
1/m or rs are close to the Horizon scale. Note that we 
choose 1/m and to be physical distances, i.e. we do 
not allow them to be comoving distances. 

B. Phenomenology 

Modifications to the Poisson equation as described 
above translate into an alteration of gravity above a char- 
acteristic scale (1/m or Ts). To illustrate this we choose 
to focus on the linear perturbation theory and compute 
the linear growth rate of structure. We do not attempt at 
building a fully covariant theory of gravity (a task well 
beyond the scope of this paper) and therefore we con- 
sider that the evolution of the background metric will 
be the one of the ACDM model currently favored by the 
data (See however [63,[6l] for a discussion of background 
evolutions). Our motivation to do so relies on the fact 
that any alternative theory of gravity will have to re- 
produce the now well observationally established ACDM 
background evolution. It seems thus fair to keep using it 
as an effective description of the evolution of the back- 
ground metric. However, if our test of gravity on 10 Mpc 
scales or so were to detect any deviation, a more acute 
theoretical description would be required. 

A modification to the Poisson equation leads to a 
change on the linear growth rate of perturbations. If 
D is defined so that S{]s.,t) = £)(k, t)6(k, tp) where to is 
an arbitrary original time (sgI . [ssj . then it satisfies the 
equation 

D + 2HD = ^^^^^f{k)D (10) 

where " ' " denotes time derivative and where we restrict 
ourselves to time independent / and k is the comoving 
wave number. The Hubble parameter, H, and its evolu- 
tion are defined as for the ACDM standard case by 



H\a) = H^(^^ + n^^ . (12) 

Solving Eq. [TO] numerically is straightforward. To fix 
boundary conditions, by analogy with the standard case, 
we consider the limit a — > and look for a solution 
D = a". Using the fact that H'^ —>■ HIVLo/o^ (we tend to 
Einstein-de Sitter) when a — > we find that the growing 
mode solution is given by 

n = i ( + 24/(fc) - l) if 1 + 24/(fc) > 0. (13) 



Numerical solutions to this equation at z = for vari- 
ous parameters are plotted in Fig.[Tl or more precisely the 
growth rate for modified gravity divided by the expected 
growth rate of structures for a ACDM model. Let's dis- 
cuss the Yukawa model first (left panel). As visible in 
Eq. m the Yukawa potential introduces a damping on 
gravity for r ^ 1/m which entails a slow down of struc- 
ture growth. More generally, an increase or a decrease 
of gravity is expected, depending on the sign of a. This 
is visible in the left panel of Fig. [1] for various as and 
ms. By construction, our models converge toward stan- 
dard gravity at the smallest scales (/ ^ 1 when k — > oo) 
and tends to a weaker but scale independent gravity for 
— <C 1 hence the flat limit for low k. The amplitude of 
this plateau depends only on a. The on-set of the tran- 
sition is however controlled by both a and m and, for a 
given a, the smaller m the lower in k the transition oc- 
curs. For a given m however, the smaller |a|, the lower in 
k is the departure from standard gravity. Those results 
agree with previous published results, as e.g. Stabenau 
and Jain 14511 and the analytical estimates of Sereno and 
Peacock |44| . 

The second model we consider is the UB model, as 
defined in Eq. [9l It only depends on the cut-off scale 
Tg and as expected, the growth of structures for scales 
larger than will be slower as compared to ACDM as 
the gravity gets weaker. This is clearly visible in the 
second panel of Fig. [TJ 



III. ANALYSIS METHODOLOGY 
A. Observables 

Making use of the previously derived linear theory pre- 
dictions, we now constrain our models using two probes 
of the growth rate of large scale structures: cosmic shear 
and galaxy clustering. 

For the former we will use the latest Canada-France- 
Hawaii- Telescope Legacy Survey (CFHTLS [ill) observa- 
tions as in Fu et al. 66] , a work extending the analyses of 
the previous release [67, 68, 69). For the latter we will use 
the Sloan Digital Sky Survey (SDSS [8^) matter power 
spectrum estimated from the clustering of Luminous Red 
Galaxies as in Tegmark et al. [T^I ■ 

Since weak gravitational lensing provides a mean to di- 
rectly image the total mass distribution as a function of 
redshift, it is potentially a powerful way to constrain ac- 
curately the growth of structures and thus any theory af- 
fecting it. The data we use here are based on the recently 
published analysis of the CFHTLS- Wide survey (release 
T0003) that spreads over 57 square degrees (34.2 after 
masking). The depth of the weak lensing catalog reaches 
a magnitude of i'j^g — 24.5 and corresponds to a galaxy 
density of about 13.3 gal./arcmin^. It comprises an effec- 
tive sample of about 1.7x10^ galaxies whose shape cor- 
relation properties have been analysed by Fu et al. [6^ . 
Although five bands will eventually be available, only the 
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FIG. 1: Ratio of linear growth factor of modified gravity versus ACDM at z = 0. From left to right, the Yukawa model and 
the UB model. 



i' band was used for this analysis. The area used to pro- 
duce this data-set is about twice larger than the previous 
release of the CFHTLS (but only 35% of the total size of 
the survey) but explores much larger scales. The width 
of this area turns out to be crucial to our analysis since 
the cosmic shear is now measured from 1 arcmin up to 4 
degree, well into the linear regime. Since we are looking 
for variation in the shape of the power spectrum, this 
wide range of scales is particularly valuable. Although it 
would be highly beneficial to our project, no tomography 
measurements has been carried out with these data yet. 
Systematic effects are constrained to be smaller than sta- 
tistical errors. "B-modes" in particular are negligible at 
all scales we use. To constrain gravity, we choose two 
lensing statistics whose properties are slightly different : 
the shear E/B correlation functions, ^e, and the com- 
pensated filter known as aperture mass statistic, Af, 



with the redshift distribution given by 



17111. The former is defined as 



ap 



^E{e)^^ r dk kPMMke) 



(14) 



whereas the latter is defined as 



MUe) = ^^, j^'^^PMJKkO) (15) 



where 
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D{k,w)Ps{£/fK{t)). 

(16) 



Ps is the density power spectrum a,t z — Q and the ge- 
ometry factor W is defined as 

W{uj)^ dwn{w) — . (17 
Jw Jk{w') 



n{z) 



exp 



(18) 



and where fx is the comoving angular diameter distance. 
Both statistics involve a different weighting of the con- 
vergence power spectrum depending on whether a wide 
or narrow kernel is favored in real or Fourier space [t^. 
As illustrated in [6^, the use of two different statistics 
provide an extra consistency check for the measurements 
and their interpretation. 

We assumed here that the estimators have been prop- 
erly calibrated and that ^+ = ( in the notations of 
[il,[69[). As Fu et al. we will discard the four small- 
est angular scales when using the M^^ statistic due to 
excessive E/B mixing. We will use the exponential red- 
shift distribution defined in Eq. (TB] as Benjamin et al. 
[sgj . It contains three free parameters z^, a and (3 that 
are determined using the photometric redshift calibrated 
on the VIRMOS VLT Deep Survey [zl. The uncertain- 
ties in a and /3 are suHiciently small that we can fix a 
and (3 to respectively 0.838 and 3.43. However we still 
need to marginalize over the uncertainties, as it will 
be described in the next sub-section. The choice of the 
Benjamin et al. [6^ redshift distribution instead of the 
Fu et al. [6^ seems somewhat inconsistent and may look 
like a useless complication. It is primarily motivated by 
the huge computational gain in keeping a and (3 con- 
stant. Although the Benjamin et al. [69] distribution is 
sligthly less accurate than the power-law one used in Fu 
et al. [66], the few percent difference is not relevant for 
our purpose but it would be if we were trying e.g. to con- 
strain the overall amplitude instead of marginalizing over 
it. 

The large-scale real-space power spectrum Ps{k) is 
measured using a sample of 2x10® galaxies from the 
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Sloan Digital Sky Survey, covering about 2x10^ square 
degrees with median redshift z ~ 0.35. We follow the 
methodology described in Tegmark tt al. (70l | and use 
the likelihood code available at [13]. Note that these 
measurements of the real-space matter power spectrum 
P{k) are up to an unknown overall multiplicative bias 
factor, b, over which we will marginalize as described be- 
low. This bias is assumed to be linear, i.e. constant, in 
the range probed by the SDSS and CFHTLS-Wide data 
sets. A strong scale dependance, although unlikely on 
those scales, could potentially be degenerate with the ef- 
fects we are looking at. 

One subtlety arises when computing the density power 
spectrum required to derive theoretical expectations for 
these observables. In both cases, to define an initial 
power spectrum we use the fitting formula of Novosyad- 
lyj et al. [74| as used by Tegmark et al. [t^. We nor- 
malize it using Ag, the amplitude of density fluctuations 
at A: = 0.05 Mpc~i (see [75]). We then modify the 
growth rate as a function of k and a according to the 
numerical results of Eq. [TOl Since the galaxy clustering 
data we are using are in the mildly non-linear regime 
{k < 0.2ft,/Mpc), we do not apply any non-linearity cor- 
rections besides the redshift space distortion one as in 
Tegmark et al. [t^I. The situation is different for cos- 
mic shear. Schematically, the smallest angular scale of 
interest to our lensing data is 1 . arcmin. Neglecting pro- 
jection effects, this angular scale corresponds to about 
0.6 Mpc/h at the median redshift of our lens population 
{z — 0.5), that is around k ~ 2h/Mpc. At these scales, 
the non-linear corrections to the power spectrum are of 
order a few and thus non-negligible. However, in the 
standard ACDM model, these corrections are computed 
using fitting formulae calibrated on numerical simula- 
tions [t^, |73| that use standard gravity. Their relevance 
to the alternative theory of gravity we consider is there- 
fore not obvious. This loop hole has already been tackled 
in the literature and does not turn out to be a critical 
issue. Stabenau and Jain [i^ used numerical simulations 
with modified gravity and found that in the context of 
our Yukawa model, both prescriptions give reasonable 
fits. Peacock and Dodds [7y| seems to provide a better 
fit for negative a (our conventions), whereas for positive 
a, both prescriptions are as good. These results have re- 
cently been confirmed and extended to include the Uzan 
and Bernardeau [55[ model we consider [47[. As a pre- 
cautionary measure, all our results will use both prescrip- 
tions. 



B. Likelihood methodology and priors 

We set up to constrain these alternative theories of 
gravity using maximum likelihood methods. We choose 
to maximize the likelihood in a five (four) dimensional 
parameter space. Two (one) of these parameters are 
the prime physical parameters we are interested in, be 
it either {a,m) in the Yukawa case ((rs) in the UB 



case). We then choose an amplitude parameter Ag cor- 
responding to the normalization of density perturbations 
at /c=0.05Mpc^^. This parameter is tightly constrained 
by WMAP 3 year measurements only, but using WMAP 
as a prior seems inappropriate in our case. The CMB 
constraints come from much larger scales than the up- 
per limits of CFHTLS-Wide and SDSS data, typically a 
few hundreds of Mpc/h. Since our modified gravity mod- 
els affect gravity and so the matter power spectrum on 
scales larger than a characteristic scale - be it 1/m or 
Vg- that far exceed the 0.1-10 Mpc range of the data, 
it is not suitable to apply a CMB prior on Ag. In- 
stead, it is preferable to normalize the power spectrum on 
smaller scales. Such a normalization is available through 
the joint measurement of the Lya forest flux fluctuations 
with other constraints on the linear matter density power 
spectrum e. 5. SDSS [78]. Seljak et al. [79| showed that 
at the 2a level, the WMAP 3 years and small scale Lya 
constraints are consistent. As such, we will use conser- 
vatively a prior on Ag derived from the WMAP and Ly^ 
joint measurements : we assume a uniform prior for Ag in 
between the 99.7%CL of WMAP 3 year only Q : 2.988 < 
log(10i°As(A: = 0.002 Mpc"^))< 3.324. Note that be- 
cause of the altered growth of structures, the scaling as 
with Ag will depend on (a, m) or (r^). This choice of nor- 
malization does not matter for the SDSS data, because 
of the unknown bias, but is more critical for weak- lensing 
observations. 

Finally, we have to take into account two extra nui- 
sance parameters. For weak lensing observables, we 
marginalize over the parameter Zg of the source distri- 
bution (Eq. lisp assuming a uniform prior within the 2a 
constraints of Benjamin et al. 69] : Zg — 1.172 ± 0.026. 
We also marginalize the SDSS likelihood over the bias 
parameter assuming a uniform prior 0.25 < b < A which 
comprises most reasonable solutions (56j . 

The limited range in scales of our observations leads us 
to restrict ourselves to the following range for our prime 
parameters. We explore —0.5 < a < 0.7 and 0.001 < 
TO < 30 in the case of Yukawa, and 0.001 < < 10 in 
the case of UB. 

Finally, we fix the other cosmological parameters to 
their best fit flat ACDM WMAP value only (Spergel 
et al. [H). Namely we consider h = 0.732, uj„i = 0.1277, 
uJb — 0.0223 and Ug = 0.958. We did not explore the de- 
pendance of our results with regards to those parameters 
although we expect it to be small given the small uncer- 
tainties with which those parameters are now measured 
(53 ■ Those parameters essentially fix the background 
evolution and the initial power spectrum. 

We explore the likelihood in this five dimensional space 
using a regular linear gridding of this hyper-volume (ex- 
cept for the TO dimension that we explore using a logarith- 
mic binning), with typically 100 samples per direction. 
Our code makes use of the existence of fast dimensions 
b that require very little computation to be explored so 
that we can explore this full hyper- volume in 10 hours or 
so using 8 1.6 GHz Intel Xeon cpus. 
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FIG. 2: Likelihood contours at 68% and 95% confidence levels for the 1/m and a parameters of the Yukawa type modification to 
gravity. The left panel corresponds to the CFHTLS-Wide constrains while the right panel corresponds to SDSS LRGs. Colored 
contours correspond for CFTHSL-Wide to the use of the M^p statistic with the halofit non-linear prescription. The dashed 
lines were obtained using the M^p statistic with the Peacock and Dodds [t^ prescription whereas the dot-dashed lines were 
obtained with the halofit non-linear prescription but using the statistic. The agreement between these various prescriptions 
and statistics is a satisfying of robustness of our measurement. As expected given the wider area covered by SDSS (42 times 
bigger than the current status of CFHTLS-Wide), the SDSS constraints are much narrower despite the bias uncertainty. 



IV. RESULTS AND DISCUSSIONS 



Fig. [2] displays the 2D 68% and 95% contour lev- 



els for the Yukawa model 



and a parameters us- 



ing CFHTLS-Wide (left panel) and SDSS LRGs (right 
panel). Other parameters have been marginalized over. 
Focusing on the CFHTLS panel first, we display several 
set of contours. The colored contours correspond to the 
M^p statistic (Eq.[in]) with the Smith et al. [tJ halofit 
prescription to model non-linearities. The long dashed 
lines corres pon d to the M^^ statistic with the Peacock 
and Dodds [70] prescription. The dot-dashed lines corre- 
spond to the statistic (Eq. [H]) with the Smith et al. 
|77| halofit prescription. To the level of accuracy we 
are interested in this work, this figure shows that those 
various methods and measurements are consistent with 
one another. This is a reassuring statement regarding 
the robustness of our constraints. We will from now on 
quote numbers from the M^^ statistic using the halofit 
prescription (colored contours). 

In the case of Yukawa type models, for a given 1/m the 
deviation from gravity increases with increasing \a\ (see 
Fig. [T]) hence it is expected constraints be centered on our 
favored value a = 0. Therefore, Fig. [2] first shows that no 
deviation from standard gravity is favored by CFHTLS 
data. Since weak-lensing constrains small scales best 
(small 1/m), we observe narrower constraints at lower 
1/m. The broadening of the contours at high 1/m shows 
the limits of our data in terms of large scale sensitivity. 
Below this broadening, say below 1/m < 0.5 Mpc, we can 
see that a and 1 /m are decorrelated, i. e. the constraints 
on a seem independent of 1/m and vice versa. This can 



be understood as follows. For a given 1/m, the transition 
to modified gravity is set and a will modify the amplitude 
of the effect on the linear growth rate as visible in the 
right panel of Fig. [TJ Since for weak-lensing the relevant 
quantity is a weighted projection of the linearly evolved 
power spectrum (with non-linear corrections applied for 
each z considered), we expect a degeneracy between a 
and the overall amplitude of the power spectrum set by 
Ag. This is illustrated in Fig. [3] where we plot the 2D con- 
tours in the a — Ag plane for 1/m < O.lMpc. The range 
of As shown here corresponds to our uniform prior on As 
motivated in section IIII Bl Clearly this choice of prior 
is a key for weak-lensing measurement and tightening it 
more would also strengthen our constraints on a but we 
choose a fairly conservative prior in order to secure the 
consistency with the latest CMB and Li/a observations. 

The right panel of Fig. [2] shows the same constraints 
using the SDSS LRGs. Clearly, they are much tighter. 
The reason for this difference can be understood in the 
following way. Since we marginalize over the amplitude 
of the power spectrum within our As priors (and also on 
the bias in the case of SDSS), our test of gravity basically 
relies on constraining the shape of the power spectrum. 
Because the SDSS data we are using cover about 2400 
square degrees with accurate photo-z for all of those (and 
a median redshift of z ~ 0.1), the power spectrum we are 
using covers the range 0.02 h/Mpc< k <0.3ft./Mpc. On 
the other hand, CFHTLS-Wide data, makes uses of only 
57 square degrees so that the errors on the same mode 
are at least 6 times larger due to cosmic variance only. 
In addition, since cosmic shear is sensitive to a weighted 
projection of the matter power spectrum (see Eq. ll6p . the 
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FIG. 4: ID likelihood distribution for a when one sets 1/m = 
0.04Mpc. Likelihood corresponds to either CFHTLS-Wide or 
SDSS. Vertical lines correspond to 68% and 95% confidence 
levels. We can see that for this scale, CFHTLS-Wide con- 
straints are only ~2 times worse than SDSS ones. This can 
be understood from the fact that lensing is sensitive to the 
overall normalization of the power spectrum. 



FIG. 3: Likelihood contours at 68% and 95% confidence levels 
for the normalisation parameter As and a parameter of the 
Yukawa type modification to gravity. We consider for this 
plot only 1/m < O.lMpc so that we explore the lower part 
of the left plot in Fig. [2] This plot illustrate the degeneracy 
between a and As that is expected given Fig. [T] where we 
can see that for a given m and redshift, a will change the 
amplitude of P{k) above a given scale. 



shape information is partially erased due to projection 
effects. This effect could be greatly alleviated if we were 
able to separate the lenses into various z planes, i.e. to 
perform tomography, but such a work is still in progress 
for the CFHTLS survey. 

On the other hand, for a given scale m well probed 
by both cosmic shear and galaxy clustering data, say 
TO =25 Mpc"^ (1/to =0.04 Mpc) CFHTLS performs bet- 
ter than one would expect from simple cosmic variance 
arguments as is visible in Fig.[31 This figure shows the ID 
marginalized likelihood distribution for a when one sets 
1/to =0.04 Mpc. It can be understood as a horizontal 
slice in Fig. [2] for 1/to =0.04 Mpc. We see in this figure 
that the 68% uncertainty level for a using CFHTLS-Wide 
is only 2.2 bigger than the SDSS ones. This comes from 
the fact that lensing does not suffer from the bias un- 
certainty so that within our amplitude prior, lensing is 
sensitive to the overall amplitude whereas galaxy clus- 
tering is not. On the other hand, as illustrated above, 
if we broaden our normalization prior then the lensing 
constraints will broaden too. 

Fig. [5] displays the ID 68% and 95% confidence levels 
on from the UB model obtained either with the cosmic 
shear or galaxy clustering observations. Once again, we 
do not see any evidence for any deviation from standard 
gravity on the scales we probed. Not surprisingly, we 
find that as for the Yukawa model, SDSS performs bet- 
ter than CFHTLS in constraining by about one order 
of magnitude. Once again, this stems directly from the 
wider range of scales probed by SDSS. 

Our results for the Yukawa model are consistent with 



those previously published in the literature using the 
SDSS results |56|. ISTl or a simpler form of cosmic shear 
measurements [54| . Our constraints on the UB model are 
in agreement with those derived using the shear 3-point 
function (3^ . 

Before discussing further our results, we reassess the 
theoretical hypothesis underlying our work. First we as- 
sume that the background evolution is identical to the 
ACDM one, which might over-estimate the effects on 
large scales [531 ■ Second, we assume that the relation be- 
tween the deflection angle and the gravitational potential 
that underlies the convergence power spectrum definition 
in Eq. [11] is unchanged as compared to GR. Models were 
this relation does not hold are as such neglected in our 
approach [H, [H, [H, |4^. Third, we consider deviation 
from the usual cosmological Poisson equation only and 
restricted ourselves to some parametric form. Those lim- 
itations come from the fact that we choose not to work 
within the frame of a well defined theory of gravity. As 
such, our approach may be described as a test of the 
Poisson equation in a ACDM cosmology. 

We did not study explicitly the dependence of our 
constraints as a function of other parameters. Sealfon 
et al. [svj showed that their variations within currently 
observationally allowed bonds had minor effects in this 
context. Besides, the current constraints are such that 
most of the major degeneracies of the ACDM are bro- 
ken. The likelihoods are thus almost Gaussian (see Fig. 
10 of Spergel et al. [l|) and marginalizing over a param- 
eter is very close to setting it at its most likely value. In 
contrast, however, the effects of massive neutrino would 
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FIG. 5: Ts likelihood distribution for the UB model. Vertical 
lines correspond to 68% and 95% confidence levels. Blue lines 
correspond to CFHLS-Wide constraints whereas red lines cor- 
respond to SDSS. The blue dashed line corresponds to the 
use of the Peacock and Dodds J7(| prescription to model non- 
linearities while the solid line corresponds to the Smith et al. 
[TTI halofit prescription. 



be degenerated with gravity ones we explore since they 
also modify the slope of the power spectrum at high k 
[43 |. On the other hand, as explained before, we also did 
not quote any constraints in terms of erg since this num- 
ber is obviously very model dependent in our case and 
a comparison with the standard ACDM value would not 
be meaningful. 

The main conclusion of our work is that we did not 
find any sign of deviation from gravity using large scales 
probes of the matter perturbations on cosmological scales 
(between 0.04 and 10 Mpc) at low redshift. Whereas 
probes of gravity on those scales at low redshifts are 
scarce, there is another important one we did not use, 
namely the Integrated Sachs- Wolfe (ISW) effect 
However, since our modifications of gravity are totally 
ad-hoc and do not correspond to any covariant theory, it 
is unclear yet how meaningful would be our extension to 
the scales relevant for ISW and we leave that for future 
work. 

Another indirect implication of our work is the fact 
that the dark matter clustering evolution as probed by 
cosmic shear measurements, i.e. the deflection of lights 
and the clustering of luminous matter as probed by 
galaxy surveys provide consistent constraints on theories 
of modified gravity on 0.04-10 Mpc scales. This agree- 
ment has some more quantitative applications that we 
leave for future work. Note also that within our models, 
our constraints imply that galactic scale physics is un- 
changed. Whereas the Yukawa model we discussed was 
originally introduced as an alternative to dark matter on 
galactic scales, the extrapolation of our constraint sug- 



gests that galactic physics is still ruled by GR. 

It is worth noticing that in the context of the lin- 
earized theory of gravity, the addition of a Yukawa term 
corresponds to setting a non-zero mass to the graviton 
fH, [6ll, [8l|. Our current results do not allow us to 
place constraint on 1/m if we leave a vary freely. If 
we were to set a = 1 as in Choudhury et al. [sl], then 
we obtain the following constraints on the mass of the 
graviton 1/m > 5.71Mpc (95% CL) using CFHTLS, and 
1/m > 67.6Mpc using SDSS; that is m < 1.11 x IQ-^^eV 
using CFHTLS and m < 9.40 x lO'^^eV using SDSS. 

Further improvements on tests of gravity on cosmolog- 
ical scales should explore more physical assumptions as 
well as estimators, like tomography. SDSS data show the 
sky coverage is of primary importance, while CFHTLS- 
Wide data reveal this survey still covers too narrow a 
part of the sky and should go one order of magnitude in 
terms of angular scales. This will be possible soon, with 
the advent of next generation large scale cosmic shear 
surveys, such as e.g.DES, DUNE, LSST or SNAP [s^. 
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